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NONLINEAR MEMBRANE SOLUTIONS FOR SYMMETRICALLY LOADED 


DEEP MEMBRANES OF REVOLUTION 

By John N. Rossettos 
Langley Research Center 


SUMMARY 


An explicit closed form solution is given for the equations of a nonlinear 
theory for deep membranes of revolution. In the theory, it is assumed that 
strains are small and rotations are moderately small, that is, of the order of 
the square root of the strains. The solution is valid when the nonlinear 
Behavior is confined to Boundary- layer regions near edges. General symmetric 
surface loads can Be conveniently included in the form of the solution. 

Results are obtained for the particular cases of a pressurized membrane and a 
rotating spherical membrane, Both attached to hubs.. Also, results are given 
for a cap of arbitrary shape attached to a cylindrical pressure vessel and for 
the arrest of a moving mass by a pressurized spherical membrane. 


INTRODUCTION 


In the design of space vehicles, membranelike structures (i.e., structures 
with negligible bending stiffness) continue to fine wide application. Large 
displacements of such structures must be considered, for instance, in the design 
of deployable space power configurations, radar antennas, and pressure vessels. 
Therefore, the solution of membrane shell problems, especially in the nonlinear 
deflection range, is a useful tool for establishing rational design procedures. 

Several derivations of systems of equations appropriate for the analysis 
of a membrane of revolution in the nonlinear range are available. (See refs. 1, 
2, and 3 .) The analysis in the present paper is based on an approximate system 
of equations given by Sanders in reference 3> in which the strains are assumed 
to be small and the rotations are moderately small. 

For satisfying the common boundary conditions Imposed in practice, it is 
known that linear membrane theory is inadequate. In such cases a linear bending 
theory solution is often added to remedy this difficulty. This situation gives 
rise to the usual edge zone where bending action insures satisfaction of the 
necessary boundary conditions. However, for membranes or very thin shells of 
negligible bending stiffness this approach is not realistic, but a nonlinear 
membrane theory can yield reasonably accurate results and avoids the greater 
complications of nonlinear bending theory. The order of the system of equations 



of nonlinear membrane theory, in contrast to linear membrane theory, is found 
to be high enough so that the usual boundary conditions can be satisfied. This 
result was first pointed out by Bromberg and Stoker in reference k, who also 
noticed that nonlinear membrane behavior is confined to a narrow region near 
the boundary. 

Nonlinear membrane theory has already been used successfully to solve par- 
ticular membrane problems. Ey using an appropriate nonlinear theory, some deep 
membrane problems with specific shapes and loading conditions have been con- 
sidered. (For example, see refs. 2 and 5 «) Also, a linearized large deflection 
approach in the case of shallow membranes has been given in reference 6 in which 
shallow spherical and conical membranes internally pressurized and fixed at 
edges are considered. 

The purpose of the present paper is to obtain a general closed form solu- 
tion to the equations of nonlinear membrane theory. The solution has not 
appeared before and is valid for deep membranes of revolution of arbitrary 
shape under symmetric surface loads. For displacement boundary conditions, the 
constants of integration are given by closed form expressions. For stress 
boundary conditions, the constants of integration can be calculated in a 
straightforward manner as indicated in the particular problems treated. 

Complete solutions are found for four particular problems. The first two 
involve a pressurized spherical membrane and a rotating spherical membrane, 
both attached to hubs. The other solutions concern a cap of arbitrary shape 
attached to a cylindrical pressure vessel and the arrest of a moving mass by a 
pressurized spherical membrane. The solutions obtained by the approach used 
are shown to be accurate to within small errors of a stated order of magnitude. 


SYMBOLS 


a characteristic dimension of membrane 

constants of integration 


E 


f 


1 > 


g 

q,S 

H 

h 


Young's modulus 


., f, 


'1 


/functions of £ which appear in general solutions 


J 


horizontal component of edge load 
thickness of membrane 
energy parameter 
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M 

N l> N e 

p 

P H ,Pv 

Pn^t 

R 

R L 

R |' R 0 

r 

s 

t 

u 

V 

V 

w 

W 1 

w l,max 

X 

z 

a 

<x T 


e i> e 9 


rigid mass 
stress resultants 


characteristic pressure loading of membrane 
horizontal and vertical components of surface load 

normal and tangential components of surface load 

radius of spherical membrane 
radius of cylinder 

principal radii of curvature 


horizontal distance 

arc length 

time 

horizontal displacement 
vertical component of edge load 
velocity 

vertical displacement 
vertical displacement at hub 

maximum vertical displacement at hub (see eq. (B35)) 


distance along cylinder 
vertical distance 

factor for arc length along meridian 
value of a at juncture of pressure vessel 
pa 

small parameter- - — 

Eh 


strains 


0 


angle in circumferential direction 



£ 

Si 

Sj 

h 

p 

cp 

♦ 

^e'^e,max 

*L 

a 


Poisson 1 s ratio 

curvilinear coordinate along meridian 

value of £ at each edge of membrane; i = 1,2 

value of £ at juncture of pressure vessel 

lower limit on certain integrals which appear in solutions 

mass density of membrane material 

angle between axis of revolution and normal to middle surface 

rotation of the normal to the middle surface in a meridional 
plane 

boundary-layer rotation, and maximum value 
rotation calculated from linear theory 
frequency of rotation 


Subscripts: 

e boundary layer 

J juncture 

max maximum 

A prime over a symbol denotes differentiation with respect to 

A bar over a symbol denotes a physical quantity. When a bar does not 
appear, the quantity is nondimen sional. (See eqs. (17).) This rule does not 
apply to the physical quantities a, E, h, and p. 

FUNDAMENTAL EQUATIONS 


The geometry of the undeformed middle surface of a general membrane of 
revolution is shown in figure 1 . Also shown is the notation for the membrane 
stress resultants along the meridional and circumferential directions N^ 

and Nq, respectively. The positive directions of the displacements u and w 

and the rotation \|r, the rotationally symmetric loads per unit area p^ and Py 

( or ^n Pf)* e(ige stress resultants H and V are given In the 

figure. 


k 



The middle surface is defined by the parametric equations 


? = r( | ) 


z = z(|) 


( 1 ) 


The parameter £ is the curvilinear coordinate in the meridional direction. 
Arc length on the undeformed middle surface is given by 


ds 2 = a%£ 2 + r 2 d0 2 


( 2 ) 


where 

1/2 

r ’ = a cos 


and primes denote differentiation with respect 
dinates are used, and the principal curvatures 

l _ _l_ 

h ~ a 

The governing equations are based on a nonlinear theory for thin shells 
given by Sanders in reference 3 in which small strains and moderately small 
rotations are assumed. The system of equations of the nonlinear membrane theory 
are derived by utilizing the principle of virtual work with appropriate strain- 
displacement relations given in reference 3 . In this manner equilibrium equa- 
tions with a consistent set of boundary conditions can be derived. 

In terms of the notation used, the strain-displacement relations are 


9 


z 1 = a sin cp 


to £ . Lines of curvature coor- 
are given by 


sin cp 
r 


(3) 


i. = —(u 1 cos cp + w f sin cp) + ~ 
£ a v / 2 


G 0 ” ? 


cos cp - u' sin cp) 
a' ' 


The stress -strain relations are 


*1 - s( H s • vNe ) 


w 

(5) 

(6) 


(7) 
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( 8 ) 


f 9 ' £(% ' *' N l) 

The principle of virtual work (which equates the internal virtual change 
in strain energy to the external virtual work of the loads pg and py and 

the edge stresses H and V) is expressed in the following form: 




+ N e 



ae ae 



Sw]a? d9 dg 


+ 



+ 



'*2 _ 
r d9 

lh 


(9) 


If the relations in equations (4) to (8) are now used in equation (9), and the 
indicated variations are carried out, integration by parts gives 



cp) - Ilj&u + JjNg^sin cp + \|r cos cp) - vjsw^ ^ 


cp - sin cp) I - aNg + arpg 1 5u d9 d| 


r d9 


>)]' 

Oj + a?Pyj 


V 


cp + t cos cp j + arpy > 5w d6 d£ - 0 


( 10 ) 


By equating to zero the coefficients of the variations of the displacements in 
equation (10), the following equilibrium equations and boundary conditions are 
obtained: 


Equilibrium equations: 



cp + cos 



-arp v 



cp - t sin 



aN 0 = -arp H 


( 11 ) 


( 12 ) 
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Boundary conditions: 

•p* 

u 


w 


At an edge or prescribe 

or N|(cos 9 - ? sin 9) = H 
or N|(sin 9 + ip cos 9) = V 


( 13 ) 

m 


The surface loading may "be divided into two general classes. In one type 
of loading, the directions of the load components py and p^ are always 

parallel and perpendicular, respectively, to the axis of revolution of the mem- 
brane; the magnitude does not change with deformation. For instance, in the 
case of a membrane spinning about the axis, i>y = 0, and p^ is given by the 

centrifugal forces. In the second general loading situation, the loads are 
fixed in magnitude but are always normal p Q and tangential p^ to the middle 

surface during the deformation. In this case, for the moderately small rota- 
tions considered, the components perpendicular and parallel to the axis of 
revolution are given by 


Pg = p t (cos cp - ^ sin 9) + P n (sin cp + t cos cp) 
Py = p^sin cp + if cos 9) - p n (cos 9 - i|r sin 9) 


( 15 ) 

(16) 


For example, for a pressurized membrane, p^ = 0 and p n equals the internal 

pressure. In the present theoiy, it is noted further that the terms with \|r 
in equations (15) and (l6) can be dropped and still be within the accuracy of 
the present solution. 

It is convenient to introduce the following nondimensional variables and 
parameters where bars denote physical quantities: 


( K 5' N 8> V ' H )- E | S (i' ? ,N 9 ,V,h) 


(w*) ’ §( 5 t’V*) 

(r,a,R|,R e ) = i(r,a,R ? ,R e ) 
(P H ,Py,P t , p n ) = i(p H ,P v ,P t ,P n ) 






(IT) 
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where e = is a small parameter which is of the order of magnitude of the 

strains. The quantity p is seme characteristic loading which in a given case 
shall he selected so that it has units of force per unit area. The quantity a 
is some characteristic dimension of the membrane. 

In terms of these nondimensional variables the governing equations are 


JrN|(sin cp + ei|r cos 

9 jj = -o-rPy 

(18) 

JrN|(cos 9 - ei|c sin 9) 

. 1 

- ccNg = -arpg 

(19) 

e, = i(u' cos 9 + w' 

b <X 

sin 9) + i ei|r 2 

(20) 

€ q = - 
0 r 


(21) 

i|r = i(w* cos 9 - 

u' sin 9) 

(22) 

e | = N | - vN 0 

e e = N e - v h 

(23) 


From equations (20), (21), and (22) the following compatibility equation may be 
obtained: 


( re eV i p 

cos cp - t sin 9 = * — ■ ■ + — e^^cos 9 

This equation will be used in the next section to obtain a solution for the 
rotation i|r. 


(210 


General Solution of the Equations 


A closed form solution of the governing equations (eqs. (l8) to (23)), 

which is accurate to within errors of order of magnitude o(e^/ 2 ) with respect 
to unity, is obtained by using boundary -layer methods and asymptotic integration. 
(See, for example, refs. 8, and 9*) At the start, assumptions which involve 
orders of magnitude are made which are verified later by the solution. Near an 
edge where boundary or continuity conditions are prescribed, it is assumed that 
a boundary layer exists. Within the boundary layer it is further assumed that 


the rotation 


t = o(e -1 / 2 ). 


and 


= o(e _1 / 2 )i = C^e -1 ), J 


i|f d| = 0(1), and 
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y = 0 (l). The stresses Ng and are assumed to be 0 ( 1 ) everywhere. 

Outside the boundary layer, = 0 (l). 

The solution is to be obtained first by solving the equilibrium equations 
(eqs . (l 8 ) and (19)) for N| and N Q in terms of the rotation \|f, and then by 

substituting these results via equations (23) into the compatibility equation 
(eq. ( 24 )) to obtain a single equation for f. Equation (l 8 ) is integrated to 
give 


rN|(sin cp + e\|f cos cp) = c-g (0 


( 25 ) 


where c is a constant of integration and 


g(S) 


-r 


arpy d| 


( 26 ) 


Since 


e\|r is small, equation (25) is used to write N| as 


Ne = C _ e|r cot cp) + 0 (e) (27) 

5 r sm cp 

This step limits this method to values of cp that are not too small (i.e., 
those for deep membranes) so that cot cp does not became too large. It is 
also noted that the term with et|r in equation (27) is retained because it will 
be differentiated later, and therefore its order of magnitude will be increased. 

Now, in order to calculate Nq, it is necessary to differentiate rN^; 
this term is given by 


( 



= -(c - g)j 


cos cp cos cp N 


sm 


2 cp 


sm 


3p 


arpy 

— — — -f- 

sm cp 



(28) 


The term Nq is obtained by using equations (19), (27), and (28) and is 
written as 

N 0 = ~ ~q S V f~ + cot T + rPfl + C^e 1 / 2 ) (29) 

a sin^tp J 

The stress-strain equations (eqs. (23)) are used next to write the compatibility 
equation (eq. ( 24 )) in the form 
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^N| - vNg^cos cp - \(f sin cp = ^ J* ) ^rN| ^ + i ei|r 2 i 


cos cp 


(30) 


If the expressions for the stress resultants are now substituted into equa- 
tion ( 30 ), the following equation for the rotation i)r is obtained: 


(c - g)cos cp (c - g)cos cp /q 
r sin <P a sin^p 


(I + *') + 


rp v cos cp cot cp - rpg cos cp 


(rp s - rp v cat cp)' - «**) - 



- vr^py sin cp + p^ cos cp^ = i|r sin cp + i ei|r 2 cos cp + o(e^/ 2 ) 
The solution for i|r in equation (jl) can be written as 


(31) 


♦ = + e + + L 


(32) 


The term ^ is the result of linear membrane theory and is of the order of 
magnitude unity- The solution for can be obtained readily from equa- 

tions (l8) to (23) by setting e = 0- The first term in equation (32), 
namely ty e , is of the order of magnitude and represents a homogeneous 

edge solution which decays rapidly away from an edge- 

When equation (32) is substituted into equation (31) various terms of 
linear membrane theory are canceled identically by the linear expression 
for ^i-e., is subtracted out at this stage j so that what remains is 

given by 


<1^. 


O 

a sin cp 


+ \ ei|r e 2 cos cp 


f c - g > 
Rt sin^cp 


■of ^e' + °( el/2 ) 


(33) 


The assumption is now made that the membrane geometry is sufficiently smooth so 
that differentiation does not increase the order of magnitude of the quantities 
a., r, and R^. Furthermore, if previous boundary-layer assumptions are 

invoked and only large terms of order o(e - -*-/^) are retained in equation (33) 
(which has a relative error of o(e^/ 2 ) compared with unity), the following 
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boundary- layer equation for \|r e 


is obtained: 


V + e = 0 


(3*0 


where 


g^sin^cp 
r(c - g) 

The solution of equation (34) can be found by noting first that it is a differ- 
ential equation which contains a large parameter, namely, l/e. Such equations 
are discussed in reference 8. Accordingly, it is convenient to make the fol- 
lowing transformations for dependent and independent variables. Let 



P = qWe 


dx = 


(35) 


Equation (34) becomes 


* 2 P 


- [j + Pl (x)]3 = o 


(36) 


where 


p ( X ) = is" 

p l l ' 4 J 16 q 5 

The quantity p-, (x) is continuous and order of magnitude unity l this holds for 

d 2 q \ ' 

smooth loading so that q" = — ^ = 0(l) ) and can therefore he neglected with 

d£ 2 1 

respect to l/e. This approximation is well within the accuracy of the present 
solution. The solution for p, then, is clearly of the exponential type and 
taking the solution, which decays away from the edge gives the following 

houndary-layer solution for \(r e : 


= dq -1 A(|)exp^- j= 


f 1 T^COd! 


(37) 


The assumption that \|r g is o(e -1 /^J implies that the constant of integra- 
tion d in equation (37) is o(e“-*-/^). 
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Now if the membrane has two edges it is assumed that they are sufficiently 
far apart so that the boundary-layer solutions can be treated independently of 
each other- (Note that, in general, d can take on a different value for the 
boundary layer near the other edge ) Once the rotation is calculated with 
equation (32), all other quantities that are expressed in terms of \Jr can be 
determined. In particular the expressions for the stresses N^ and Nq are 

given by equations (27) and (29). 


The displacements u and w are now to be determined. For brevity the 
symbol (i = 1 or 2) will represent the value for either of the edges 

or By using equations ( 21 ) and (23) together with the expressions for the 

stress resultants N| and Nq, the displacement u near edge can be 

determined: 


u 


rfi(l) 

q. 


- rS(|) + 


/id 

~q$F 


rf 2^) ex P 





+ o(e 1 / 2 ) (38) 


where 


^(1) 


q^sin cp va^sin^tp 

rR g + 


f 2 U) = ^L£lLiP 


S(£) = rpy cot cp - rpjj 


The first two terms in equation (38) represent the linear membrane solution 
for u, and the third term (a decaying exponential) is a boundary-layer contri- 
bution due to the rotation \[f e . 


From equations ( 20 ) and ( 22 ) and from previous relations the following 
expression for w* is obtained: 



Equation ( 39 ) can integrated to give 
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w = cf ^(0 - f^(|) + J at cos cp d£ + o(€ 1 / 2 ) ( 40 ) 

where 

* i) -l iI( ?+ ^N dl 

fl,u> ' K I + ' “ vS 8in cp ) ai 

The constant of integration in equation (40) has been incorporated into the 
lower limit of the integrals which appear. If the expression for the 

rotation f (given by eqs. ( 32 ) and (37)) is used in equation (^0), the dis- 
placement w near the edge can be written as 


v = cf 5 (t) + f“ 5 (i) — 

Q. 

where 

*•> ■ r a cos cpt L d| - f^Cl) 

% 

fg(t) = a cos cp 



The first two terms in equation (4l) are the result of linear membrane theory. 

In arriving at equation (4l), the facts that the boundary layer is exponentially 
decaying and that q(|) and fg(i) are smooth functions (i.e., q'(i) and 

fg'(i) are 0(l)) have been used to evaluate approximately (within the accuracy 

of the present solution) the following integral: 



f 6 U)* e 41 



-M) 


(te) 
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with 


The lower limit |jj of the integral in equation (42) is replaced by oo 
negligible error because of the rapid decay of the integrand and because 

has a value which is outside the boundary layer (given in the section on bound- 
ary conditions ) . 


Boundary Conditions 

If the boundary conditions are prescribed with respect to displacements, 
and if the membrane is attached to the hubs at each edge £ ^ (i = 1 , 2 ), the 

displacements must vanish there. If the two hubs are equal in size and the 
membrane shape and loading are symmetrical about a midpoint between the hubs, 
then the lower limit in the integrals associated with equation (4l) is 

taken to be the value of £ at this midpoint. From symmetry, it is clear that 

w (^n) = °* 

The constants of integration c and d are now evaluated explicitly for 
the class of problems where boundary conditions are prescribed on displacements. 
If the membrane of revolution is attached to hubs, the boundary conditions at a 
hub, say at £ = £^ ^the hub at £2 can "^ e treated similarly and independently 

of £^, are given by 


N - 0 

(43) 

N - 0 

(44) 


The solutions for the constants of integration are given for the case of sym- 
metry between hubs since they will be used for the particular problems treated 
in appendix A. The extension to the case of no symmetry between hubs, although 
not considered in this report, offers no new difficulty. In the symmetrical 
case it is sufficient to consider only one hub at £ = £-^. The two boundary 

conditions (eqs. (43) and (44)) yield two equations for the constants c 

and d. From equations ( 38 ) and (4l) the following two conditions are obtained: 


and 


fl( g l) ^ df 2 ( g l) 
a ( E l’°) + q3/ 4 ( Sl , c) 



of 3 ( g l) + ?5(!l) 


/ed 

<l 3/k (h’ c ) 



= 0 


(45) 


(46) 
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When equations (45) and (46) are solved, c and d are given hy 


c 


(sf 6 - f 5 f 2 )r f 2 2 sin <p - gf x f 6 
rf^fg^sin cp - f x f 6 


-4=tl 


( 47 ) 


1 

q3A(cf 5 + 

/i 

L f 6 J 


(48) 


Note that the order of magnitude of d is for smooth loading. Equa- 
tion (37) then shows that the assumed order of magnitude of is confirmed. 

With c and d thus determined by equations (47) and (48) all relevant quan- 
tities may be calculated. Stresses and displacements are calculated in appen- 
dix A for particular problems by specializing the foregoing general results. 


If the boundary conditions are prescribed with respect to stresses, the 
following equations must apply at an edge: 


N|(cos cp - ei(r sin cp) = H 

(49) 

N|(sin cp + e\|r cos cp) = V 

(50) 


These equations are the nondimens ional form of equations (13) and (l4). The 
stress resultants given by H and V cannot be prescribed independently of 
each other but must be related by the following expression which is obtained 
from equations (49) and (50): 


V cos cp - H sin cp = 0 


) 


(51) 


Equation (51) is necessary in order to insure that rotations are not greater in 
magnitude than the order of the square root of the strains, since this is the 
basis for the present theory. However, for the physical problems involving 
deep membranes which actually occur in practice, the behavior implied by equa- 
tion (51) is the usual case. This is clarified by the cases treated in appen- 
dix B, which presents more specific indications of actual stress boundary 
conditions . 
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RESULTS AND DISCUSSION 


Particular problems are solved in the appendixes which illustrate the use 
of the general results obtained. Stresses and displacements are calculated in 
appendix A for the special cases of a spherical membrane under internal pres- 
sure and of a rotating spherical membrane, both attached to hubs. These calcu- 
lations are performed by specializing the general formulas in the text. The 
results of these calculations, which give stresses and displacements, are shown 
in figures 2 to k. 

In figure 2 the maximum physical rotation \jr e max at the hub is plotted 

against hub angle g n for several values of, the small parameter e = — . In 

Eh 

the parameter e, R is the radius of the membrane, whereas the quantity p 
represents the internal pressure for the pressurized membrane and the centrif- 
ugal force per unit area (p = phO%) for the rotating membrane. In the case of 
the pressurized membrane, t e ,max is independent of the hub angle g 1 whereas 
in the case of the rotating membrane it decreases with Figures 3 a/nfl 14 - 

show the nondimensional stress resultants Ng and N g and the displace me nt w 
in each case, which occur near a hub having a hub angle g^ = k^° . The differ- 
ences between linear and nonlinear membrane theory are clearly indicated by the 
boundary- layer regions near the hub. The curves for the displacement w show 
that linear theory does not allow satisfaction of the boundary condition at the 
hub ^i.e., w(g-^) = 0 ^ so that nonlinear membrane boundary layers are needed 

here to accomplish this condition. 

Although the stresses are positive for the pressurized spherical membrane, 
it is interesting to note that, in the case of the rotating spheric al membrane, 
the circumferential stress N g becomes negative in a small region near the hub 

so that wrinkling can occur in that region. 

In appendix B, two problems are given for which boundary conditions are 
prescribed on both stresses and displacements. In the first problem, juncture 
stresses are calculated for a cylindrical pressure vessel. Equations (52) 

(eqs. (B19) in appendix) give the result of computing the membrane stress 
resultants Ng and N 0 at the juncture of an arbitrarily shaped cap and a 

cylinder. The stresses are 




N e “ 1 " 55 




(52) 


The results in equations (52) show that Ng is always positive, whereas N g 

can become negative when the ratio of the meridional radius of curvature of the 
cap at the juncture to the radius of the cylinder is less than l/k ^i.e., 

Rg j < It is interesting to note that the result of linear membrane theory, 
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obtained by dropping ef e ' in equation (55) (eq. (B5) in the appendix), indi- 
cates that Nq becomes negative when ^ j is only less than l/2. 


N Q = - 




2a 


sin%> 




') + 


sin cp 


(55) 


Also, in contrast to nonlinear membrane theory, linear theory yields discontin- 
uous stresses at the juncture. These stresses are shown in figures 5 6 for 
the case of a spherical and a torispherical cap, respectively. The stress Nq 

is plotted in these figures, with the boundary layers necessary to maintain con- 
tinuous stresses. 


The second problem in appendix B concerns the arrest of a moving rigid 
mass attached to a pressurized spherical membrane. This problem shows that the 
general results in the text offer a convenient means of obtaining the necessary 
relation between the unknown edge stress resultant and deflection at points of 
contact between the mass and membrane (i.e., relations similar to influence 
coefficients). Maximum stresses and displacements are calculated for a range 
of values of initial velocity of mass. The results are shown in figure 7 where 
the quantities are plotted against a kinetic energy parameter 


k E = 


Mv 2 / 1 

€pr 5^« sin 


«1 


) 


CONCLUDING REMARKS 


Closed form analytical solutions have been obtained for the equations of a 
nonlinear membrane theory in the case of deep membranes of revolution under 
arbitrary but otherwise axially symmetric surface loads where the membrane 
geometry and loading are also assumed to be smoothly varying. The solutions 
are found to be accurate to within errors of the order of the square root of e 
with respect to unity, where € is a small parameter which has a magnitude 
characteristic of strains. The nonlinear behavior is represented by exponential 
type terms which decay rapidly away from the boundary. 

Expressions for the constants of integration which appear in the solutions 
have been obtained in general form for the class of problems in which boundary 
conditions are prescribed on displacements. The general results can be special- 
ized in a straightforward manner to handle specific problems. This specializa- 
tion is accomplished in the particular cases of a pressurized membrane and a 
rotating spherical membrane, both attached to rigid hubs. The stresses and 
displacements for these two examples exhibit the expected boundary-layer behav- 
ior. Also, whereas the stresses are always positive in the pressurized mem- 
brane, the circumferential stress can become negative in the rotating membrane 
for certain values of the parameter e, so that wrinkling can occur. 
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For problems in which boundary conditions are prescribed on stresses, it 
is found that, according to the present theory, the vertical and horizontal 
edge stress resultants cannot be prescribed independently of each other but 
must be related in a particular way. Two examples are given in which boundary 
conditions are prescribed on both stresses and displacements. The first con- 
cerns juncture stresses for a cylindrical pressure vessel for which boundary 
layers eliminate the usual stress discontinuities characteristic of linear mem- 
brane theory. Also, in this case, nonlinear theory predicts that negative cir- 
cumferential stresses at the juncture can be avoided if the ratio of the merid- 
ional radius of curvature of the end cap to the radius of the cylinder at the 
juncture is greater than l/4 . (Linear membrane theory gives 1/2 for this 
ratio.) The second example involves the arrest of a moving rigid mass attached 
to a pressurized spherical membrane. For this example, the general solution 
provides a convenient linearized form for the necessary relation between edge 
stress and corresponding edge deflection (i.e., similar to influence coeffi- 
cients). This allows a straightforward determination of maximum stresses and 
deflections for a range of initial velocities of the mass. 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Administration, 
Va., November 22, 


1965 . 
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APPENDIX A 


PARTICULAR PROBLEMS WITH BOUNDARY CONDITIONS ON DISPLACEMENTS 


Two examples are given for deep memLrane problems in which boundary condi- 
tions are prescribed on displacements only (i.e., the membranes are attached to 
rigid hubs). The solutions are obtained by a straightforward specialization of 
the general results given in the text. 


Pressurized Spherical Membrane Attached to Hubs 

In the special case of a spherical membrane attached to a hub and of rad- 
ius R and internal pressure p, the various quantities defined for the general 
problem are specialized. From equations (l), (2), ( 3 ), and (17) and from other 
relations the quantities are reduced as follows: 


R|=Rq = 1 a = 1 cp = £ 


a = R 


P n = P = 1 


r = sin | 


p. = 0 



From equations (15), (l6), and (17) 


Pjj = sin £ + ©|r cos | 


(A2) 


Py = -(cos | - a|/ sin |) 


(A3) 


The solution for the rotation given by linear membrane theory ^ can be found 
from equations (l8) to (23) by setting e = 0. In this case 


= 0 


(A4) 


The various other relevant quantities which appear in the general solution given 
in the text may now be evaluated. 
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g(l) = - | sin 2 | 
fiU) = 1 + v 

f 2 d) = 1 

f,d) = f 1 . + e d| = (1 + v)log tan -| 

f^d) = - f - - V sin | d| = 1 - 2 - V cos 5 

J rt/2 

f 5 (S) - di-^l COS | 
fg(l) = COS 5 

S(|) = 1 


"N 


\ (A5) 


The constants of integration c and d given by equations (47) and (48) can 
now be calculated by using equations (A5). The solution for c is 


c = 


, fl + V , ,1 - v\ 

cos Ht" - 1 + -g-) 

(l + v)^log tan - cos csc 2 £-^ 


= 0 


(A6) 


Then from the definition of q(|) obtained from equation (34), q(|) = 2 so 

that the solution for d is 


d = 


1_ 1 - v 

/e ( 2 ) x /4 


The maximum rotation at the hub £ 


li is (see eq. (37)) 




e,max 



1 - v 

/2i 


(A7) 


(A8) 
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The solutions for the stresses and displacements are determined from the gen- 
eral results given by equations ( 27 ), (29), (37)> (38), and (^l). For this 
problem they are 



N e 


1 _ i_^ 

2 2 


t = 


-d - v) 

yie 



u = 



w 



cos £ + 



cos 







J 


(A 9) 


The solution given by equations (A9) is accurate to within errors of order 

o(e^/^) with respect to unity, as was indicated by the general results in the 
body of the paper. Also note that the linear membrane solution is exactly what 
remains when the exponential terms are dropped from equations (A9). 


Rotating Spherical Membrane Attached to Hubs 

For a rotating spherical membrane the radius of the membrane attached to 
a hub is also denoted by R. However, the characteristic loading quantity p, 
which is discussed in relation to equations ( 17 ), arises from centrifugal 
forces and is given as 


p = phH R 


(A10) 


where p is the mass density of membrane material and H is the frequency of 
rotation. The relevant quantities used for the rotating spherical membrane are 
given as 


“ R 0 “ 1 

a = 1 

-0 

li 

irr* 

a = R 

p = phfl 2 R 

% = 


r = sin | 


e 


pQ 2 R 2 


E 



(All) 
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The linear solution for the rotation \|r L is obtained from equations (l8) 
to (23) by setting e = 0. The result is 

t L = -(3 + v)sin £ cos | (A12) 

Under these conditions the functions of | are given by 


g(l) = 0 

f x (|) = 1 + v 

^(i) = 1 


f 3 (l) 


= f -r + — > ■ = (l + v)log tan i 

J n / 2 sin | 2 


y v sin5| d| = -v cos | + v cos^l 

rt /2 


f\(l) 

f 3 (|) = cos5| + V COS I 


= COS I 

S(|) = -sin 2 i 

The constant c given by equation (47) then reduces to 

cos |q_ 

c = 

o £l 

cos csc^^ - log tan 

Por this case, from equation (34) , 


L (^l,c) = 


sin 2 ^ 


> (A13) 


(Al4) 


This result, together with equation (48), gives the following result for d: 
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(A15) 


d = 


The expression for tJt is now 

g y max 

c(l + v)csc | 3 _ - sin^^ 

^e,max 

The stresses and displacements are then given by 



(A16) 


N t = 


1 sin 2 | 


Nn = - 


® sin 2 | 


+ = .max ex P 


c ( ri l cos £ 1 " cos M\ . o k 

- / - sin l x .in t * e ,ma* exp 1 — Lj + sin I 


lcos - cos | 

- — - (3 + v)sin | cos £ 

-fee 


> (AIT) 


— c ( 1 + v) / ec sin *1 , 
U = — i + 1/ : If 


sin | 


sin £ T e,max 


exp 


cos - cos £ 

L + 


-/i5 


sm- 


w = c(l + v)log tan — 4- cos3| + v cos £ 



The form of these equations and the numerical results based on these equations 
are essentially the same as those given in reference 2. 
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APPENDIX B 


PARTICULAR PROBLEMS WITH BOUNDARY CONDITIONS 
ON BOTH STRESSES AND DISPLACEMENTS 


End Cap Attached to a Cylinder Under Internal Pressure 

Stresses and displacements are calculated at the juncture of an end cap of 
arbitrary but rotationally symmetric shape and a cylinder under internal pres- 
sure. Solutions are first obtained separately for the cap and cylinder by 
using the general results found previously. Then the constants of integration 
are evaluated by using suitable matching conditions at the juncture. For this 

problem a boundary layer is expected at the juncture ^i.e., cp = fig. 5 j • 

Cap solution . - The internal pressure p is taken as the characteristic 
loading quantity of equations (17). For the cap solution then, p n = 1 and 

p^ = 0 so that the horizontal and vertical loads are 


P H = sin cp + a|r cos cp (Bl) 

Py = -(cos cp - ety sin cp) (B2) 

Also, for the cap solution, the quantity a used in equations (17) is taken as 

the radius of the cylinder R^ so that the small parameter is given by 

= Eh; 

6 Eh 


Row, from equation (26) and with the relation r' = a cos 9 

r 2 

g(S) = - \ 


(B3) 


Equation (27) is next used to obtain an expression for . The constant of 

integration c is taken to equal zero so that the solution reduces to the 
linear membrane solution outside of the boundary layer. The stress result- 
ant N| is then given by 


N. = 

5 2 sin cp 


(B 4 ) 


2k 



From equation (29) Nq is obtained: 


N e = - 


2a sin' 


+ *•’) + ^ 


(B 5 ) 


The boundary- layer solution for the rotation \|r e is given by equation (37) and 
in this case becomes 


where 


t e = dq _1 /^(|)exp 



q^UJdl 


(B6) 


l(l) = 


2a^sin5p 

73 


and | T is the value of | at the juncture . The displacements can be obtained 

O 

with the help of equations (38) and (4l) and are given by 


u 




sxn ep \ 



W = fy(l) - 





(B7) 

(B8) 


where 


f T (?) ■ Iq (“ 008 ’*1 + IT + 2B*"ln" - avr ) 45 (B9) 

JC 

Cylinder solution .- For the cylinder solution, a = 1 , r = 1 , and cp = — . 

Note also, for this solution, that the coordinate | is the ratio of the dis- 
tance along the cylinder to the radius R^. The corresponding solutions for 

stresses and displacements are obtained and, to an accuracy with an error of 
o(e 1 /^) compared to unity, give 


N | = c 2 


(BIO) 
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N 0 = 1 - c 2 eV 

(Bll) 

u = 1 - C 2 ^ e T - vc 2 

(B12) 

w = (c 2 - v^i 

(B13) 

= exp (- 7^ l e - *jI) 

(B14) 

The constants of integration are now evaluated by requiring that the displace- 
ments are continuous and the forces are in equilibrium at the juncture of the 
cap and cylinder. The general form for the boundary conditions is given by 
equations ( 13 ) and (1*0, so that the matching conditions are as follows. At 

the juncture, where cp = these conditions require that 


u cyl = ^ap 

(B15a) 

w cyl = w cap 

(B15b) 

(* H l)= y l “ (* H l)cap 

(Bl6a) 

Ns = N. 

s cyl 5 cap 

(Bl6b) 


The subscripts cyl and cap refer to the cylinder and cap. The constant c^_ 

appears only in equation (B15b), so that the other three conditions can be 
solved directly for d, c 2 , and d^. Equations (Bl6a) and (Bl6b) yield 


1 

e 0 = — 
2 2 






j 


(B17) 


Equation (Bl5a) is then used to get 


d = 


1 



(B18) 


26 



where R| j and aj denote values of R^ and a at the juncture . If the 
stresses are now computed at the cap-cylinder juncture, it is found that 
and N 0 are given by 


N 





(B19) 


Arrest of a Moving Mass by a Pressurized Spherical Membrane 

The fact that masses may be attached to pressurized spherical membranelike 
structures during deployment of space vehicle packages makes it of interest to 
calculate the maximum stresses which would occur in the arrest of such masses. 
The general solutions in this paper offer a convenient means for calculating 
these stresses. In the present problem it is assumed that, in the deployment 
operation, no extensional strains are suffered by the folded membranelike 
material as it unfolds until the spherical shape is completely formed; it is 
also assumed that the velocity of the moving mass is known (from other consid- 
erations) at this instant. From then on membrane strains exist and contribute 
to the arrest of the moving mass. 

For the pressurized spherical membrane (see fig. 8), 


= R 0 = 1 

a = R 

From equations 


a - 1 


p n = p 

(15), (16), 


cp = I 


p t = 0 


r = sin | 

£ = PR 
Eh 


and. (IT) 


p y = -(cos | - ajr sin £) 


p H = sin | + eif cos £ 


From equation ( 26 ) 


g(i) 


sin 2 £ 

2 



(B20) 


(B21) 

(B22) 


(B25) 


The edge of the membrane which is attached to the rigid mass M will suffer 
no displacement u while the mass is moving (fig. 8). At the same time there 
will be a vertical edge force V because of the deceleration of the mass, so 
that the appropriate boundary conditions at | = |-j_, in dimensionless form, are 


27 


I 



"(Si) - 0 

N|^sin + ©|r cos = V 
Equation (B23) and equation (2j) yield 


(B24) 

(B25) 


c + 


• 2a 

sir? 


- 


sin 2 £ 


-(l - ei|r cot |) + 0(e) 


(B26) 


If this expression for N| is substituted into equation ( B 25 ) , the terms 
O^e 1 / 2 ) cancel identically so that the following relation is obtained: 

V = — + — sin £, + 0(e) 

sin 2 1 ' 

Next, an expression for N 0 is obtained from equation (29): 


(B27) 


n 9 = 1 


b * ; la2 % - jajA exp JL r t 

sln2| l -/i J Sl 


q^dl 


(B28) 


where 


<1 = 


sin 2 £ 


c + 


sin 2 £ 


If the boundary condition (eq, (B2^)) is now applied and terms of order 

are neglected since they are small compared to unity, the following 
equation is obtained: 


u 


( 8l )--(^ + 5)[ 1 + v -^ lA (^)] = 0 


(B29) 


Equation (B29) can be used to calculate the constant d: 

iA 

i /c + i sin2|, ' 

1 + vl 2 1 


d = 




sin 2 ^ 


(B30) 
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An expression for the displacement w at I = will he needed. First, note 
that 


s(i) = -i 


fjd) = (1 + v)log 

^U) = hl-H cos | 


tan 4 
2 




(B31) 


Equation (4l) evaluated at | = gives 


w 


(il) - A( ?1 )V - B( ?1 ) 


(B32a) 


or 


v.!i + 2 

A A 


(B32b) 


where 


ii 

tan — 


B 


A (ll) = (1 + v)^sin l x log 
(ll) = i|(l - v)cos !]_ + (l + v)sin 2 i 1 log tan 


.) 


Equations (B32) provide a relation between the as yet unknown V and the dis- 
placement w at £ = Another relation is needed to determine the maximum 

deflection which occurs during the arrest of the mass. This relation 

can be obtained by equating the work done by the edge resultant V to the 
kinetic energy of the moving mass at the instant the spherical shape is formed. 
When written in terms of physical quantities this relation is 


r w l,max _ _ ^ 

-2jtR sin £, / y dw, = ± Mv 2 

x J 0 2 


(B33) 


Equation (B33) can he written in terms of nondimens ional quantities: 


-X 


w- 


l^max 


V dwi - *e 


(B3*0 
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I 



where the nondimens ional parameter kg 


is seen to he 


1 Mv 2 ' 

sin €pr3 

proportional to the initial kinetic energy of the mass. The integration in 
equation (B3^) can he performed directly if the relation between V and w^ 

given hy equations (B32) is used. The following result is then obtained: 



A (h) 



W l,max 


(B35) 


The maximum deflection w-, can he determined from equation (B35)* Equa- 

jl y max 

tions (B27), (B30), and (B32) are associated with this deflection and can he 
used to calculate c, d, and Y } and therefore the maximum stresses. Maximum 
values of Ng, Ng, w, and at 6 = (^1 = ^5°) are plotted in figure 7 

against the parameter kg. 
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Pressurized membrane 




= ^5°) for pressurized and rotating 
(v = 0.3). 
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